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ly^ . Abstract 

o, 

O . Integration of the form J^ f{x)w{x)dx, where w{x) is either sin(a;x) or cos{ujx), is 

widely encountered in many engineering and scientific applications, such as those 

(^ ■ involving Fourier or Laplace transforms. Often such integrals are approximated 

H I by a numerical integration over a finite domain (a, b), leaving a truncation error 

equal to the tail integration Jl^ f {x)w{x)dx in addition to the discretization error. 
This paper describes a very simple, perhaps the simplest, end-point correction to 
approximate the tail integration, which significantly reduces the truncation error 
and thus increases the overall accuracy of the numerical integration, with virtually 
no extra computational effort. Higher order correction terms and error estimates for 
the end-point correction formula are also derived. The effectiveness of this one-point 
correction formula is demonstrated through several examples. 
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tion error. 



1 Introduction 

Integration of the form J^ f{x)w{x)dx, where w{x) is either sin{ujx) or cos(a;a;), is 
widely encountered in many engineering and scientific apphcations, such as those involv- 
ing Fourier or Laplace transforms. Often such integrals are approximated by numerical 
integrations over a finite domain (a, b), resulting in a truncation error /^°° f{x)w{x)dx, in 
addition to the discretization error. One example is a discrete Fourier transform (DFT), 
where there is a truncation error due to cut-off in the tail, in addition to the discretization 
error. 

In theory the cut-off error can always be reduced by extending the finite domain at the 
expense of computing time. However, in many cases a sufficiently long integration domain 
covering a very long tail can be computationally expensive, such as when the integrand 
f{x) itself is a semi-infinite integration (e.g. forward Fourier or Laplace transform), 
or when the integrand decays to zero very slowly (e.g. a heavy tailed density or its 
characteristic function). Much work has been done to directly compute the tail integration 
in order to reduce the truncation error. Examples include nonlinear transformation and 
extrapolation (Wynn 1956, Alaylioglu et al 1973, Sidi 1980, 1982, 1988, Levin and Sidi 
1981) and application of special or generalized quadratures (Longman 1956, Hurwitz and 
Zweifel 1956, Bakhvalov and Vasileva 1968, Piessens 1970, Piessens and Haegemans 1973, 
Patterson 1976, Evans and Webster 1997, Evans and Chung 2007), among many others. 
This paper describes a very simple, perhaps the simplest, end-point correction to account 
for the tail integration over the entire range (6, oo). The treatment of the tail reduces 
the usual truncation error significantly to a much smaller discrete error, thus increasing 
overall accuracy of the integration, while requiring virtually no extra computing effort. 
For the same accuracy, this simple tail correction allows a much shorter finite integration 
domain than would be required otherwise, thus saving computer time while avoiding 
extra programming effort. To our knowledge this result is not known in the literature and 
we believe it deserves to be published for its elegant simplicity and broad applicability. 
Though it is possible that our formula is a rediscovery of a very old result hidden in the 
vast literature related to numerical integration. 

The paper is organized as follows. In Section 2, we derive the tail integration ap- 
proximation and its analytical error. A few examples are shown to demonstrate the 
effectiveness of the tail integration approximation in Section 3. Concluding remarks are 
given in Section 4. 



2 Tail integration 

Consider integration J^ f{x) sin{ux)dx. Without loss of generality, we assume u = 1 
(a change of variable y = ux results in the desired form). For J^ f{x) cos{ux)dx the 
derivation procedure and the resulting formula are very similar. In the following, we 
assume that 

• The integral J^ f\x) sm.{(jjx)dx exists; 

• All derivatives /'■^•'(a;) exist and — )> as A; — )■ oo. 

2.1 Piecewise linear approximation 

The truncation error of replacing j^ f(x) sm(x)dx by /„ f{x) sm(x)dx is simply the tail 
integration 

0< 

f{x) sm{x)dx. (1) 



b 

For higher accuracy, instead of increasing truncation length at the cost of computing time, 
we propose to compute the tail integration It explicitly by a very economical but effective 
simplification. Assume f{x) approaches zero as x — ?> oo and the truncation point b can be 
arbitrarily chosen in a numerical integration. Let b = Nir, where A^ is some large integer. 
Dividing integration from Nn to oo into cycles with an equal length of tt yields 

/ f{x)sm{x)dx = '^Ik, Ik = / f{x)sm{x)dx. (2) 

Nn ^=^ {N+k)7T 

Now assume that f{x) is piecewise linear within each 7r-cycle, so that each of the integrals 
Ik in (2) can be computed exactly. That is, in the range [(A^ + k)7i, {N + k + l)n], we 
assume that f{x) is approximated by 

f{x)^fk + ifk+i-fk), (3) 

n 

where fk = f{{N + k)TT). Substitute (3) into (2), then analytical integration by parts of 
each Ik in (2) gives 



/"^ oo oo 

fix) sm{x)dx = J2lk^Y: i-lf ^'{fk + fk+i) = (-1)^/0 = i-irfiNn). (4) 



This elegant result given by (4) means that we only need to evaluate the integrand f{x) at 
one single point x = Nir (the truncation point) for the entire tail integration, replacing the 
truncation error with a much smaller round-off error. As will be demonstrated later, this 
one-point formula for the potentially demanding tail integration is remarkably effective 
in reducing the truncation error caused by ignoring I^- 

2.2 Higher order correction terms and error estimates 

Formula (4) can be derived more generally through integration by parts, and a recursive 
deduction gives us higher order correction terms and thus error estimates. Integrating (1) 
by parts with b = Nn, we have 

OO OD 



f{x)sm{x)dx = {~iy^f{NTr)+ / p^>{x)cos{x)dx, (5) 

Nit Ntt 

where f^^\x) = df{x)/dx. If we assume f{x) is linear within each 7r-cycle in the tail, 
then the integration j^^ f^\x) cos{x)dx vanishes, because within each vr-cycle f^^'{x) is 
constant from the piecewise linear assumption and J^^ cos{x)dx = for any integer 
k, and f^^\oo) — )> as /(oo) — )> 0. Thus, under the piecewise linear assumption, (5) and 
(4) are identical. Continuing with integration by parts in (5) and noting f^-^\x) — )■ at 
infinity, we further obtain 

oo oo 



f{x) sm{x)dx = {-ly^fiNn) - / p'>{x) sm{x)dx, (6) 

Nw Ntt 

where /'■^•'(x) = d^f{x)/dx'^. Equation (6), as well as (5), is exact - no approximation 
is involved. The recursive pattern in (6) is evident. If we now assume that the second 
derivative /'■^•'(x) is piecewise linear in each 7r-cycle in the tail, then (6) becomes 

f{x) sm{x)dx ^ (-1)^ (/(iVvr) - f^'\Nn)) . (7) 

Nn 

With the additional correction term, (7) is more accurate than (4). In general, without 
making any approximation, from the recursive pattern of (6) we arrive at the following 
expression for the tail integral 



fix) sm{x)dx = i-lffiNn) + ^ (-1)^+7(2^) (iVvr) + (-1)^ / /(''^^(x) sm{x)dx, 

Ntt *=1 Ntt 



where k > 1, f^'^'^^Nn) is the 2k-th order derivative of f(x) at the truncation point. As 
will be shown later with examples, typically the first few terms from (8) are sufficiently 
accurate. The error in using formula (4) is readily obtained from (8) 

fc-i °° 

^T = E (-1)^+7^'^^^^^) + (-1)' / f^''\x)sm{x)dx, (9) 

In deriving (8), we have assumed all derivatives exist and f^'^\oo) = 0. Under certain 
conditions, the infinite series in (8) and (9) represents the integral asymptotically as 
A^ — 7- oo, i.e. we have the asymptotic expansion 

/ fix) sm{x)dx ~ E i-lf^"f^'\N7r) as N ^ oo 

Nn *=0 

For example, if we assume that, for some 7 < 0, 



/(™)(a;) =0(x''-™) as A^^oo, m = 0,l,2,..., 

then the integral term on the right-hand side of (8) can be bounded by CkN^'"^^^^ as 
A^ — )■ oo, for some positive constant C^, and the series converges to the integral. 

The derivatives approaching zero as a; — > oo is a consequence of the existence of 
integral (1). Otherwise, if f'^^^oo) ^ 0, integral (1) does not exist, which is evident 
form (5). Applying this argument recursively, all derivatives f^''\oo) = 0, if they exist. 
Obviously if f{x) is a power function (e.g. l/^/x), the ratio f^^\x)/ f(x) is of the order 
1/x as a; — )■ oo, so is the ratio f'^'^^^\x)/ f^^\x). This implies that, for a power-like 
function, each error term in (9) decreases by two orders of magnitude from its preceding 
term as the index number k increases by one. 

Remark. Note that there is no truncation error in (4) and the error is a discretization 
error in nature. In theory, the tail integration error can be estimated by (9). In practice, 
however, derivatives of integrand at the truncation point may only be evaluated numeri- 
cally. The assumption of piecewise linearity, although reasonable for f{x) at large x, may 
appear to be rather crude for a high precision computation. However, we recall that we 
are only trying to reduce the already small truncation error It and a reasonable approxi- 
mation in It could lead to significant improvement in the overall accuracy of integration. 
For example, suppose a relative error of 1% due to ignoring truncation and 10% error in 
evaluating the tail integration using the very simple formula (4). The overall accuracy 
with this tail integration added is now improved from 1% to 0.1% (1% times 10%). This 
improvement by an order of magnitude is achieved by simply evaluating the integrand at 



the truncation point. The assumption of a piecewise hnearity apphes to a broad range 
of functions, thus the special tail integration approximation can have a wide application. 
Note, piecewise linear assumption does not even require monotonicity - f{x) can be os- 
cillating, as long as its frequency is relatively small compared with the principal cycles in 
sin(x), as demonstrated in one of the examples below. 

If the oscillating factor is cos(a;) instead of sin(x), we can still derive a one-point 
formula similar to (4) by starting the tail integration at (A^ — l/2)7r instead of Nn. In 
this case, the tail integration is 



fix) cos{x)dx ^(-l)'V {{N - l/2)7r) . (10) 

(Af-l/2)7r 

Also, the tail integration approximation can be applied to the left tail (integrating from 
— oo to —Nn) as well, if such integration is required. 

It is known from the literature that truncation is better at extrema of the oscillatory 
part than at the zeros (Lyness 1986, Espelid and Overholt 1994 and Sauter 2000). Trun- 
cating at xt = [N — l/2)7r, the extrema for sin{x), we obtain an expression for the tail 
integration or the truncation error similar to (8) 



fc-i 



Jf{x)sm{x)dx = (-l)^/«(xr) + ^(-l)^+7(2^+i)(xr) 

CO 

+ {-l)^ f f(^'^){x)sm{x)dx. (11) 



XT 



The leading term of the truncation error is now f^^^x^) in (11), compared with f{xT + 
l/27r) in (8). Assuming \f^^\xT)/ /{xt + l/27r)| < 1 for some large Xt, e.g. when f{x) is 
a power-like function, then it is obvious truncation at extrema has a smaller truncation 
error. However, our formula is about the reduction of the truncation error by including an 
approximation of the tail integration. If truncation is done at x = (A^ — l/2)7r instead of 
X = Nil, then the first correction term will be (— l)^+^/*^^)((A^— l/2)7r), involving the first 
derivative of f{x). In many important applications the first derivative of f{x) cannot be 
evaluated accurately. For example, when inverting a characteristic function of a compound 
distribution, f{x) itself is a semi-infinite integration of an oscillatory function, which could 
only be obtained numerically. Taking finite difference of a numerically evaluated function 
will in general reduce the accuracy by an order of magnitude. So for general purposes the 
truncation is chosen at the zeros, i.e. at a; = Ntt. 



Of course, if derivative of f{x) is in closed form and can be accurately evaluated, 
truncation and correction at extrema will indeed be more accurate, with a leading error 
term of f^^^x-r). But we could also include the second derivative term for the truncation 
at zeros, with a leading error term of f^^\xT)i and so on. In general when higher order 
derivatives can be computed precisely, then one can include some higher order terms to 
reduce truncation error further and it does not matter much whether the truncation is 
done at extrema or at zeros. 

3 Examples of tail integration 

The effectiveness of the above tail integration approximation is now demonstrated in a 
few examples. Introduce the following notations 

rOD 

Ie = fix) sm(x)dx, 

Jo 

i-ttN 

I{N-k) = / f{x) sm{x)dx, 
Jo 

/^(A^Tr) = / f{x) sm{x)dx. 

JnN 

In all the following examples the exact semi-infinite integration Ie is known in closed 
form, and its truncated counterpart /(A^tt) is either known in closed form or can be 
computed accurately. For simplicity in all the examples A^ is taken to be an even number, 
i.e (—1)^ = 1. The exact tail integration /^(A^Tr) can be computed from IxiNn) = 
Ie — I{N7[). We compare I{Ntt) + /{Nit) with /(A^tt) and compare both of them with 
the exact semi-infinite integration Ie- The error reduction can be quantified by comparing 
the "magic" point value given by formula (4) with the exact tail integration It{Ntt). Also 
note that the analytic formula for the error of using (4), Et = Ie ~ [/(A^tt) + /(A^tt)], is 
given by (9). 

Example 1: f(x) = e~°^, (a > 0). 

In this example, the closed form results are 






Ie = e "'' sm(x)dx = -, I(Nn) = / e "''^ sm(x)dx 

J 1 + a^ J 





Figure 1 compares the "magic" point value /(A^vr) representing simplified tail integra- 
tion with the exact tail integration It{Nti) = Ie — I{N-k) as functions of parameter a for 



A^ = 10, i.e. the truncated lengths It = lOvr. The figure shows that a simple formula (4) 
matches the exact semi-infinite tail integration surprisingly well for the entire range of 
parameter a. Corresponding to Figure 1, the actual errors of using formula (4) are shown 
in Table 1, in comparison with the truncation errors without applying the correction term 
given by (4). Figure 2 shows the same comparison at an even shorter truncated length of 
Air. The error of using (4) is I^tI = ^^ exp(— a7rA^)/(l + a^). 

If a is large, the function f{x) = e~"^ is "short tailed" and it goes to zero very fast. 
The absolute error {etI is very small even at A^ = 4. The relative error (against the 
already very small tail integration), given by I^tI / Ie — I{2Nn) = a^, is actually large 
in this case. But this large relative error in the tail approximation does not affect the 
high accuracy of the approximation for the whole integration. What is important is the 
error of the tail integration relative to the whole integration value. Indeed, relative to 
the exact integration, the error of using (4) is {etI / \Ie\ = (^ exp(— ayrA^), which is about 
2.7 X 10^^^ at A^ = 4. The condition /*■*'-' (x) — > as A; — )■ oo is not satisfied in this case 
if a > 1. However, as discussed above, the application of formula (4) does not cause any 
problem. 

For a small value of parameter a, the truncation error will be large unless the truncated 
length is very long. For instance, with a = 0.01 the truncation error (if ignore the tail 
integration) is more than 70% at It = lOvr (A^ = 10, as the case in Figure 1), and it is more 
than 88% at It = 47r (A^ = 4, as the case in Figure 2). On the other hand, if we add the 
"magic" value from formula (4) to approximate the tail integration, the absolute error of 
the complete integration \et\ due to this approximation is less than 0.01%, and the relative 
error is I^tI = a^ = 0.01% at both It = IOtt and It = 47r. In other words, by including 
this one-point value, the accuracy of integration has dramatically improved by several 
orders of magnitude at virtually no extra cost, compared with the truncated integration. 
For the truncated integration I{Ntt) to have similar accuracy as I {Air) + /(47r), we need 
to extend the truncated length from An to SOOvr for this heavy tailed integrand. 

Example 2: f{x) = 1/y/x . 

This example has a heavier tail than the previous one. Here, we have closed form for Ie, 

but not for I or It, 

oo Nit . , , 

Ie = I — ^dx =\ -, /(A^vr) = / — ^dx. 



'X \ Z J X X 

• 



/ or It can be accurately computed by adaptive integration functions available in many 



numerical packages. Here we used IMSL function based on the modified Clensliaw-Curtis 
integration method (Clenshaw and Curtis 1960; Piessens, Doncker-Kapenga, Uberhuber 
and Kahaner 1983). 

Figure 3 compares the "exact" tail integration /^-(A^vr) = J"^^ sm{x) / y/xdx with the 
one-point value /(A^vr). Again the one-point approximation does an extremely good job. 
Even at the shortest truncation length of just 27r the one-point approximation is very 
close to the exact semi-infinite tail integration. Applying the analytical error formula (9) 
to f{x) = l/y/x, we have 

- , lx3x... x(4/c-l) 

Taking the first three leading terms we get Et ~ —0.00695 at A^ = 2 and Et ~ 
-2.392 X 10"^ at A^ = 20. The relative error I^tI //E(2A^7r) is about 1% at A^ = 2 
and it is about 0.002% at A^ = 20. Apparently, if the extra correction term f'^'^^Nn) is 
included as in (7), the error Et reduces further by an order of magnitude at A^ = 2 and by 
several orders of magnitude at A^ = 20. Corresponding to Figure 3, the actual errors of 
using formula (4) are shown in Table 2, in comparison with the truncation errors without 
applying the correction term given by (4). 

Figure 4 shows the truncated integration /(A^tt) and the truncated integration with 
the tail modification (4) added, i.e. / + /(A^tt), along with the correct value of the full 



integration Ie = v 7r/2. The contrast between results with and without the one-point tail 
approximation is striking. At the shortest truncation length of 2tt {N = 2), the relative 
error due to truncation for the truncated integration {Ie — I{Ntt))/Ie is more than 30%, 
but with the tail approximation added, the relative error {lE — I{N'K) — f{NTi))/lE reduces 
to 0.5%. At IOOtt, the largest truncation length shown in Figure 4, the relative error due 
to truncation is still more than 4%, but after the "magic" point value is added the relative 
error reduces to less than 0.5 x 10"^. 

Another interesting way to look at these comparisons, which is relevant for integrating 
heavy tailed functions, is to consider the required truncation length for the truncated 
integration to achieve the same accuracy as the one with the "magic" value added. For 
the truncated integration I{Ntt) to achieve the same accuracy of /(27r)+/(27r) (integration 
truncated at one-cycle plus the "magic point value), we need to extend the integration 
length to TTOOtt. For I{N7c) to achieve the same accuracy of /(lOOvr) + /(lOOvr), the 
integration length has to be extended to more than IO^^tt! On the other hand, if we add 
the tail approximation /(TTOOtt) to /(TTOOvr), the relative error reduces from 0.5% to less 
than 10~^^! This error reduction requires no extra computing, since /(77007r) is simply a 

9 



number given by I/v^TTOOtt. 

Example 3: f{x) = cos{ax)/x, a < 1. 

We have remarked that the piecewise hnear assumption does not require monotonicity, 
i.e. f{x) can be oscillating, as long as its frequency is relatively small compared with the 
principal cycles. For example, when the function f{x) is the characteristic function of a 
compound distribution, it oscillates with its frequency approaching zero in the long tail. 
In the current example with f{x) = cos(Q;x)/a;, there is a closed form for Jg, but not for 
/ or It, 

/■ cosfaa;) sinfa;) , vr ?/„,^ x f cosfaa;) sin(x) , 

Ie = / ^ — -^dx =-, /(2iV7r) = / ^ — ^-^dx, a < 1. 

J X 2 J X 



Figure 5 compares the "exact" tail integration Ix{Ntt) with the one-point approxima- 
tion flNn) for the case a = 0.2. Again the one-point approximation performs surprisingly 
well, despite f{x) itself is now an oscillating function, along with the principal cycles in 
sin(a;). The piecewise linearity assumption is apparently still valid for relatively mild os- 
cillating f{x). Corresponding to Figure 5, the actual errors of using formula (4) are shown 
in Table 3, in comparison with the truncation errors without applying the correction term 
given by (4). Not surprisingly, the errors are larger in comparison with those in examples 
1 and 2, due to the fact that f{x) now is itself an oscillating function. Still, Table 3 
shows the truncation error is reduced by an order of magnitude after applying the simple 
formula (4). 

Figure 6 compares the truncated integration I^Nn) against / + /{Ntt), along with 
the correct value of the full integration Ie = n/2. At truncation length Gvr, the shortest 
truncation length shown in Figures 5 and 6, the relative error l^j^l /Ie is less than 0.06% 
and it is less than 0.01% at A^ = 100. In comparison, the truncated integration without 
the end point correction has relative error of 2.7% and 0.2%, respectively for those two 
truncation lengths. Applying the analytical error formula (9) to f{x) = cos{ax)/x and 
noting sin(aa;) = and cos{ax) = 1 with a = 0.2 and x = IOOtt, we obtain 



Et^-I + ^ + T- + ^ h x = It = IOOtt, 



a^ 2\ fa* 12^2 24' 
— + ^ + ^ + —^ 

•Xi lAj / \ •X' iXi iXj 

where only the first two leading terms corresponding to the 2"^^ and 4*^* derivatives are 
included, leading to St ~ 0.0001273 + 5.07749 x 10"^ at A^ = 50 that agrees with the 
actual error. Similar to the previous example, if we include the extra correction term 
f^'^\Nn), the error reduces further by two orders of magnitude at A^ = 100. 

10 



The purpose of Example 3 is to show that the piecewise hnear approximation in 
the tail could still be valid even if there is a secondary oscillation in f{x), provided 
its frequency is not as large as the principal oscillator. If the parameter a is larger 
than one, then we can simply perform a change of variable with y = ax and integrate 
(l/y) sin{y/a) cos{y) in terms of y. Better still, for any value of a, we can make use of the 
equality cos{ax) sin(a;) = (sin(a;— Q;a;)+sin(a;+Q;a;))/2 to get rid of the secondary oscillation 
altogether before doing numerical integration. In practice, the secondary oscillation often 
has a varying frequency with a slowly decaying magnitude, such as in the case of the 
characteristic function of a compound distribution with a heavy tail. In this case it might 
be difficult to effectively apply regular numerical quadratures in the tail integration, but 
the simple one-point formula (4) might be very effective. 

All these examples show dramatic reduction in truncation errors if tail integration 
approximation (4) is employed, with virtually no extra cost. If the extra correction term 
j(2)(^2A^7r) is included, i.e. using (7) instead of (4), the error is reduced much further. 

4 Conclusions 

We have derived perhaps the simplest but efficient tail integration approximation, first 
intuitively by piecewise linear approximation, then more generally through integration by 
parts. Analytical higher-order correction terms and thus error estimates are also derived. 
The usual truncation error associated with a finite length of the truncated integration 
domain can be reduced dramatically by employing the one-point tail integration approx- 
imation, at virtually no extra computing cost, so a higher accuracy is achieved with a 
shorter truncation length. 

Under certain conditions outlined in the present study, the method can be used in many 
practical applications. For example, the authors have successfully applied the present 
method in computing heavy tailed compound distributions through inverting their char- 
acteristic functions, where the function f{x) itself is a semi- infinite numerical integration 
(Luo, Shevchenko and Donnelly 2007). 

Of course there are more elaborate methods in the literature which are superior to 
the present simple formula in terms of better accuracy and broader applicability, such 
as some of the extrapolation methods proposed by Wynn 1956 and by Sidi 1980, 1988. 
The merit of the present proposal is its simplicity and effectiveness - a single function 
evaluation for the integrand at the truncation point is all that is needed to reduce the 
truncation error, often by orders of magnitude. It can not be simpler than that. Also, in 

11 



some applications the function f{x) may not even exist in closed form, for instance when 
f{x) is the characteristic function of some compound distributions as mentioned above, 
then f{x) itself is a semi- infinite integration of a highly oscillatory function, which could 
only be obtained numerically. In such cases some of the other more sophisticated methods 
relying on a closed form of f{x) may not be readily applicable. 
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Table 1: Error of using formula (4), St-, in comparison with the truncation error It if 
formula (4) is not applied, corresponding to Figure 1 in Example 1. 



a 


St 


It 


0.001 


9.4 X 10-^ 


0.9391 


0.01 


5.3 X 10^5 


0.5334 


0.1 


1.8 X 10"^ 


0.0018 


1 


2.6 X 10-28 


2.6 X 10-28 


10 


0.0 


0.0 



Table 2: Error of using formula (4), e-p, in comparison with the truncation error It if 
formula (4) is not applied, corresponding to Figure 3 in Example 2. 



X 


St 


It 


An 


1.0 X lO-''^ 


0.2241 


lOvr 


1.1 X 10"^ 


0.1422 


207r 


1.9 X 10-^ 


0.1006 


507r 


1.9 X 10-6 


0.0637 


IOOtt 


3.4 X 10-^ 


0.0318 
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Table 3: Error of using formula (4), Et-, in comparison with the truncation error It if 
formula (4) is not applied, corresponding to Figure 5 in Example 3. 



X 


St 


It 


207r 


4.2 X 10-^ 


0.0105 


407r 


2.1 X 10^^ 


0.0053 


eovr 


1.4 X 10"^ 


0.0035 


SOtt 


1.1 X 10-^ 


0.0026 


IOOtt 


8.4 X 10-5 


0.0021 
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Figure 1: Comparison between the exact tail integration /^ /(x) sin(x)(ix and simple 
one-point approximation f{N7f) from formula (4), when f[x) = e^"^ and N = 10. 
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Figure 2: Comparison between the exact tail integration J^ f{x) sin{x)dx and simple 
one-point approximation /(A^vr) from formula (4), when f{x) = e""^ and A^ = 4. 
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• one-point approximation 
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Figure 3: Comparison between exact tail integration /^ f{x) sin{x)dx and the simple one 
point approximation (4), /{Nn), as functions of truncated length It = Nn, 4: < N < 100, 
when f{x) = l/y/x. 
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full integration 

• truncated plus one-point tail approximation 
A truncated integration 
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Figure 4: Comparison between truncated integration I(Nn) = Jq f{x) sm{x)dx and the 
truncated integration plus the one-point approximation of tail integration, I{N7r)+f{NTT), 
as functions of the truncated length It = Ntt, 4 < A^ < 100, where f{x) = \j ^fx. 
The solid line represents the exact value of the full integration without truncation error. 



Ie = /(oo) = Jn/2. 
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— exact tail integration 
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Figure 5: Comparison between exact tail integration /^ f{x) sin{x)dx and the simple one- 
point approximation (4), f{Nn), as functions of truncated length It = Nn, 6 < N < 100, 
f{x) = cos{ax)/x, a = 0.2. 
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Figure 6: Comparison between truncated integration I{Nn) = Jq f(x) sm(x)dx and the 
truncated integration plus the one-point approximation of tail integration, I{NTT)+f{NiT), 
as functions of the truncated length l^ = Nit, 6 < N < 100, where f{x) = cos{ax)/x, 
a = 0.2. The solid line represents the exact value of the full integration without truncation 
error, Ie = -^(oo) = 7r/2. 
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